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Abstract. Recent collective flow measurements including higher moment event 
anisotropy from the PHENIX experiment are presented, and the particle type, beam 
energy dependence and the relation with jet modification are discussed. The measured 
higher order event anisotropy with event plane defined at forward rapidities and the 
long range correlation with large 77 gaps are both consistent with initial geometrical 
fluctuation of the participating nuclei. In 200 GeV Au-|-Au collisions, higher order 
event anisotropy, especially simultaneous description of V2 and V3, is found to give 
an additional constraining power on initial geometrical condition and viscosity in the 
hydrodynamic calculations. V2, V3 and V4 are almost unchanged down to the lower 
colliding energy at 39 GeV in Au-|-Au. The measured two particle correlation with 
subtraction of the measured v„ parameters shows a significant effect on the shape 
and yield in the associate particle A(/) distribution with respect to the azimuthal 
direction of trigger particles. However some medium responses from jet suppression or 
jet modification seems to be observed. Direct photon V2 has been measured in 200 GeV 
Au-f-Au collisions. The measured V2 is found to be small at high px as expected from 
non-suppressed direct photon Raa — 1, which can be understood as being dominated 
by prompt photons from initial hard scattering. On the other hand, at lower pT < 4 
GeV/c it is found to be significantly larger than zero, which is comparable to other 
hadron V2, where thermal photons are observed. 



1. Higher order event anisotropy measurements 

In order to measure the higher order event anisotropy, the event planes from the higher 
moments at various rapidities are defined with various reaction plane detectors; the 
Reaction Plane detector (RXN: \r]\ = 1.0 - 2.8), the Muon Piston Calorimeter (MPC: 
|?7| = 3.1 - 3.7), the Beam-Beam Counter (BBC: |?7| = 3.1 - 3.9) and the Zero Degree 
Calorimeter (ZDC: |?7| > 6.5). FigJUshows the event plane correlations between forward 
and backward rapidities for the same moments in panels (a) and (b), for the different 
moments in panels (c) and (d). The significant positive correlations between the same 
order higher moments are visible in panel (b), which is consistent with the initial 
geometrical fluctuation [1]. A weak negative correlation between 1*^* and 3'''^ moment 



i A list of members of the PHENIX Collaboration can be found at the end of this issue. 
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Figure 1. higher order event plane correlation between same and different harmonics 
and between different rapidities, where the subscript n=l~4 and superscript X=A~F 
indicated on $^ in the figure are the order of harmonics and the index of the following 
sub-detectors, respectively; A:RXN(+), BiBBC^"), C:MPC(+), DtMPC^"), E:BBC(+) 
and F:ZDC^+/^\ where (+) or (— ) indicates the positive or the negative rapidity. The 
descriptions of the sub-detectors and their 77 acceptances are given in the text. 



event planes can be an indication of rapidity anti-symmetric V3 contribution as seen in 
panel (d), which is further confirmed by the correlation measurements with respect to 
the 1*^* moment event plane from the spectators shown in panels (e) and (f), where the 
weak anti-symmetric V3 contribution shows the same sign as vi at the same rapidity. No 
visible correlation between 2^^^^ and 3'''^ harmonics seen in panel (d), this also indicates 
that the presence of geometrical fluctuation is the dominant effect. 

The measured V2, V3 and V4 with respect to the same order event plane \E'n defined at 
forward rapidity in RXN are shown as a function of px for several centrality slices in each 
panel in Figj2] [2] and as a function of Npart for two px bins in Figj3] p], Npart dependence 
is compared with a generalized eccentricity up to the higher moments. Significant higher 
moments compared to V2 are observed for central collisions, and smaller centrality 
dependence is seen especially in V3. One should also note that this V4(\E'4) seems to 
be larger than the previously measured V4(\E'2) [3j, because of the additional fluctuation 
contribution in \E'4 angle with respect to \E'2 angle. 

The measured V2 and V3 are compared with various model calculations in FigJH 
Since the initial generalized eccentricities e„ between Glauber model and CGC-KLN 
model have opposite trends for 2'^'^ and 3"^*^ moments because of different fluctuation 
contribution with respect to the initial overlapping geometry, the shear viscosity over 
entropy ratio rj/s used in hydrodynamic calculations will be more constrained by 
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Figure 2. measured v„ parameters as a function of px for different centrality slices 
at 200 GeV Au+Au collisions 
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Figure 3. measured v„ parameters as a function of Npart (for two different px bins) 
at 200 GeV Au+Au collisions 



simultaneous fitting to the measured V2 and V3, compared with the previous V2 alone 
fitting. The existing hydrodynamic calculations seem to describe the experimental 
measurements using the Glauber type initial condition with a smaller shear viscosity 
and the other calculations with the MC-KLN type initial condition with a larger shear 
viscosity are disfavored by the data|ll |5]. 

Identified hadron V3 is shown as a function of px for tt'^ , and p+p in FigJS] 
The similar mass ordering of V3 compared with V2 is observed at relatively lower px 
region as one might expect from the radial flow picture from hydrodynamic expansion, 
while at intermediate px region, the similar baryon - meson splitting in V3 given by 
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Figure 4. comparison of V2 and V3 with, various model calculations 
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Figure 5. Identified hadron V3 as a function of pt at 200 GeV Au+Au collisions 



the number of constituent quarks could possibly be seen just like the quark number 
scaling of V2, which might be related to the partonic degrees of freedom during the V3 
evolution. Both observations at least tell us that the V3, which is originated from the 
initial geometrical fluctuation, seems to reflect the hydrodynamic evolution as triangular 
expansion. 

The higher harmonics event anisotropy v„ measurements for charged particles and 
the identified hadron V2 measurements has also been done for both 62 GeV and 39 
GeV Au+Au collisions, the observed v„ and V2 results for given pt and centrality show 
very similar results to those of 200 GeV Au+Au collisions, which would indicate that 
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the hydrodynamical properties are mostly unchanged within the coUision energy region 
39-200 GeV[6]. 

2. Implication on two-particle correlation measurements 
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Figure 6. two particle azimuthal correlation between forward and central rapidities 



FigJH] shows the two particle azimuthal correlation between forward and central 
rapidities, where the left panels are with RXN and the right panels are with MPC for 
forward rapidity. The both top and bottom panels show the same data points (black) 
and the data are fitted with the Fourier function shown as red curves up to the 2°^ order 
for top panels and up to the 4**^ order for the bottom panels. The difference between data 
points and fitted curve is shown as blue points, where the correlation is well described by 
including up to the 4*"^ order as shown in the bottom panels, while the differences shown 
in top panels show significant 3'''^ moment signal, which has previously been considered 
as the near-side ridge and away-side double-peak like structure. These fitted parameters 
can then be used to extract the higher order event anisotropy and has been confirmed 
to be consistent with the results from the event plane method. One has to keep in 
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mind that the same experimental observations (near- and away-side structures) are just 
expressed in a different way by the higher order harmonic parameters v„. 
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Figure 7. two particle azimuthal correlation within the central rapidities 



Fig J7] shows the inclusive 7-hadron correlation measured within the central rapidity 
region (|?7| < 0.35) in central Au-|-Au collisions compared to p-|-p collisions [7]. The 
left panel shows the result with the V2 BG correction, while the right panel shows 
the corrected results with the measured V2, V3 and V4 parameters. The jet associated 
yields for Au+Au collisions are determined by the absolute normalization method 
without the ZYAM assumption. Although the away-side double-peak like structure 
has been significantly reduced with including the higher moment parameters v„ in the 
BG correction, there seems to be some remaining effects from the jet suppression and/or 
modification in terms of broadening and enhancement in both near- and away-side shape 
relative to the shape in p+p collisions. 

3. Direct photon elliptic event anisotropy measurements 

The direct photon carries various different contributions; (1) prompt photon from initial 
hard scattering at high px, which is expected to give zero V2, (2) fragmentation photon 
after the parton energy loss, which should follow the hadron as positive V2, (3) jet 
conversion photon and Bremstrahlung photon, which would have negative V2, since this 
is a part of energy loss itself, (4) thermal photon, which would carry the history of the 
dynamical evolution of quark gluon plasma. The direct photon V2 has been determined 
by the measured vr*^ and inclusive photon V2 and relative direct photon yield over the total 
inclusive photon yield including various hadronic decay photons. This measurement has 
been extremely difficult because of the relatively small direct photon signal especially at 
low pt region. Using the virtual photon measurement ^ [9] , the accuracy of the relative 
direct photon yield determination over the total inclusive photon yield has been greatly 
improved. This improvement has given a significant reduction of systematic error on 
the extracted direct photon V2 value below 4 GeV/c in px in the following figures. The 
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Figure 8. 7r°, 7'"'^- and 7'^"'' V2 parameters as a function of px at 200 GeV Au+Au 
minimum bias collisions 
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Figure 9. 7r°, 7™'^- and 7'*"- V2 parameters as a function of px for three centrality 
selections at 200 GeV Au+Au collisions 

measured 7r°, 7™'^' and 7'^*''' V2 are shown in FigJS] and Figj9] as a function of px in 200 
GeV Au+Au collisions for minimum bias and for three centrality bins, respectively|10]. 
The data are shown with two event planes defined at two different rapidities. The 
difference between two event planes is clearly seen for 7r° and 7*"'^' V2 at higher px 
region expected from the possible non-flow contributions, however the difference in the 
^dir. smaller than the systematic error bands. 

At the higher px region, the direct photon V2 is consistent with zero as also seen 
in FigJTOl which confirms the previous direct photon Raa measurement to be about 1 
because of the dominant prompt photon production. On the other hand, at the lower 
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Pt region as shown in FigJH] and FigJHl the direct photon V2 is found to be significantly 
larger than zero and similar to the hadron V2. Various model calculations assuming 
entire history of thermal photon emissions which includes a significant early time (small 
V2) contribution have failed to explain the experimental measurement. 
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Figure 10. 7r°, 7'"'^- and 7'^"'' V2 parameters (integrated above 6 GeV/c) as a function 
of Npart at 200 GeV Au+Au collisions 



4. Conclusion 

The recent results on the collective flow measurements from the PHENIX experiment 
are presented. The higher order event anisotropy has been measured to be consistent 
with initial geometrical fluctuation followed by hydrodynamic expansion. The initial 
condition and shear viscosity in hydrodynamic calculations are found to be more 
constrained by these measurements. This also strongly influences on the understanding 
of the extracted shape of two particle correlation from the jet-medium interaction. The 
direct photon elliptic event anisotropy V2 has been observed to be small as expected 
from the prompt photon production at high px, while it has been observed to be as 
large as the hadronic V2 at lower pt, where thermal photons are seen. 
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